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Abstract 

The observed limiting fragmentation of charged particle distribu- 
tions in heavy ion collisions is difficult to explain as it does not apply to 
the proton spectrum itself. On the other hand, string percolation pro- 
vides a mechanism to regenerate fast particles, eventually compensat- 
ing the rapidity shift (energy loss) of the nucleons. However a delicate 
energy-momentum compensation is required, and in our framework 
we see no reason for limiting fragmentation to be exact. A predic- 
tion, based on percolation arguments, is given for the charged particle 
density in the full rapidity interval at LHC energy (\/s = 5500 GeV). 



Recentely, the phenomenon of limiting fragmentation, [1], or longitudinal 
scaling, was rediscovered in the framework of high-energy heavy ion collisions 
[2]. In general the inclusive particle distribution, dn/dy, is a function of the 
central rapidity y and of the centre of mass energy y/s, or of A = y — j/j and 
yb, where y\, is the beam rapidity, with y/s = mi,e yb . Limiting fragmentation 
essentially means that as A becomes larger than some y^ dependent threshold 
A , A > Ao(|/6), dn/dy becomes a function only of A, 

p(A,y b ) -> /(A), (1) 

dy A>A (y b ) 
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independent of y b . As A (y b ) decreases with y b the region in A of limiting 
fragmentation increases with the energy. 

It has been argued that limiting fragmentation reflects the fast parton 
distribution in the beam [3] . However, as a sizeable fraction of the fast parti- 
cles building up the limiting fragmentation behaviour are nucleons (protons), 
[4], one requires specific parton correlations to generate the observed nucle- 
ons. This is done, for instance, in the Dual Parton Model, [5], by introducing 
valence diquarks which, directly or indirectly, [6], produce baryons, thus pre- 
serving the flow of baryon number. 

It is important to remark that the proton spectrum at high rapidity does 
not scale in the sense of limiting fragmentation, [4], presenting a shift (A) B = 
(y) b — Ub increasing in absolute value with energy, at least for intermediate 
energies [7]. Theoretically, in QCD evolution models, a transfer of energy 
and momentum from fast partons to the sea is also expected, [8]. 

In such circumstances how is it possible to obtain overall limiting frag- 
mentation? A possible solution is given by string percolation models as 
percolation implies not only a summation in colour, [9,10], but, as well, a 
summation in momentum, [11,12]. In a sense, percolation is a mechanism for 
reacceleration of particles. 

In string percolation models strings are produced along the collision axis. 
In the impact parameter plane if the area of interaction is 7iR 2 , the projected 
discs from the strings have an area nr 2 and there are N s strings, the transverse 
density parameter rj, 

v = (^m , (2) 

is the relevant parameter in percolation. If 77 1 the strings are independent 
and do not overlap, if n > 1 the strings fuse and percolate. 

When strings overlap, due to random colour summation, the particle 
rapidity density dn/dy is not the sum of the particle density fi\ of each 
string, and the average transverse momentum < pr > is not the single string 
average transverse momentum pi. In general we have 

dn/dy = F{n)N sni , (3) 

and 

< Pt >= Pi/\[F{n) , (4) 
where F(r)) is the colour summation reduction factor: 

1 1 - e -' r i 
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Note that rj -> 0, F(r}) -> 1, and that as 77 -> oo, Ffa) -> l/?? 1 / 2 . 

Regarding the momentum summation, if we have a cluster of A" strings, 
each string made up of partons with Feynman x- variable —x\ and +xi, re- 
spectively, such that yi = yb + haxi, [12], we have for the end (forward 
rapidity of the iV-cluster string 



With percolation just one cluster is formed and A" becomes the number N s 
of strings: 



A simple model can be imagined where, at relatively low energy, there is a 
transfer of momentum from the valence string (leading proton) with creation 
of identical sea strings, followed, at very higher energy, by the mechamism 
of percolation with regeneration of fast strings, [12]. One obtains a flat 
distribution in rapidity ending at y = y^ a (see (7)). 

The consequences of percolation become very clear: a decrease of particle 
density at mid rapidity, see (3), and an extension of the length of the forward 
sea rapidity distribution from ~ y\ to ~ In iV s , see (7). 

If we impose in the model energy conservation in AA collision when there 
are N part participating nucleons, 



Vn = Vi + In N . 



(6) 



Vn 3 = yi + \nN s . 



(7) 





where E(y) 



(iut) coshy and ^Js = m b e yb we obtain, as y^ — > 00, 





where 



A = In N s - y b . 
If we write for the asymptotic behaviour of the number of strings 



(10) 



N, ~ s x ~ e 



(11) 



we obtain 



A = - 



(12) 
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with 

a = l-2X. (13) 

One should notice that 1 + A is, in this approach, the intercept of the bare 
Pomeron, [12]. As const. < N s < ^Js one also sees that < a < 1. 

In order to make practical use of (9), we need to estimate the asymptotical 
behaviour of the ratio (m T )/m b . We take two simple models: 

Model I - One assumes that the relation (4) for (p T ) applies as well to 
(tut), [12]. This means that large pt physics is dominating everywhere. 

Model II - One assumes that the ratio (mr) /mj is energy independent. 
Such approximation is probably more realistic, as one is mostly considering 
the forward rapidity contributions where the px behaviour is expected to be 
almost energy independent [13]. 

Note that in (m T ) it is also implicit an averaging over rapidity. For a 
discussion of this problem, see [14]. 

We can now estimate the asymptotic behaviour of the integrand of (9), 

^-FMNsfi! — > e- Ao , (14) 

where A = — <%6 with, for central collisions, i.e., F(rj) — > 1/n 1 ^ 2 , 

a = 3/2A (Model I), (15) 

and 

a = A (Model II). (16) 
It is clear, see Fig. la, that conservation of energy requires 

-A + A = (17) 

or a = a. One further obtains 

A = 2/7 (Model I), (18) 

A = 1/3 (Model II). (19) 

In order to arrive at the asymptotic behaviour of dn/dy we simply have 
to divide the integrand of (9) by (mrj/mj, (see Fig. lb) to obtain 

dn/dy ~ e~ A o ~ e a ' Vb , (20) 
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Figure 1: a) The integrand of (9) as a function of e A . b) The normalized 
dn/dA distribution as a function of A. In both cases A = — cn/&. 

with 

a' = a - 1/2A = 2/7 (Model I), (21) 

and 

a' = a = l/3 (Model II). (22) 

The model (with two variations) that we have considered corresponds to 
a step function distribution: dn/dA = e _A <>, A < A ; dn/dA = 0, A > A 
(see Fig. lb). It trivially satisfies limiting fragmentation, (1), with /(A) = 0. 

We shall now generalize the model by introducing the Fermi distribution, 

e" A o 

2/N part dn/dA = , (23) 

with Aq = —a'yb , A = (1 — 2\)yb = —O-Ub and S being a parameter. In the 
limit 5^0 we, of course, recover the step function. In what follows we shall 
not distinguish between pseudo-rapidity and rapidity. 

In the Table we present the values of a', a and 6 in the case of Model I 
and Model II. We show as well the results for a', a and 5 when fitting RHIC 
data for N 2 dn/dy, with Eq. (23). Our results are very similar to the fits of 
[15]. We note that the experimental values for a' and a are not very different, 
as in Model II. This means that the energy dependence of (mr) /uib is weak. 

In the Table we have also included the values of a', a and S resulting from 
an overall fit to RHIC/PHOBOS data. 
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Figure 2: Overall fits to PHOBOS/RHIC data using Eq. (23) - see the Table. 
The dotted-dashed line is a prediction for LHC, ^Js = 5500GeV. 

In Fig. 2 we show our curves from the overall fit, in comparison with 
data, and our prediction for LHC, 5500 GeV data. At mid rapidity we 
expect dn/dy to be about 1500. 

We come now back to the question of limiting fragmentation. It is clear, 
from Fig. 2, that we do not have strict limiting fragmentation in the A>0 
region. 

From (23) one sees that, for A 3> A , limiting fragmentation simply 
means 

(24) 

which requires, for 5 > 0, the limiting fragmentation condition (see (23)) 

-A' + ^ = , (25) 

or 

5 = Ao/A'q = a/a' . (26) 
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The fact that we do not have limiting fragmentation is not our choice. The 
RHIC/PHOBOS data, fitted with the parameterization (23), clearly shows 
that relation (26) is not obeyed: a/a'>l and 5 ~ 2/3 < 1 (see the Table). 

In our estimates, we have assumed that at present energies and for large 
number of participating nucleons, F(rf) — > l/i] 1 ^ 2 ,!] being large enough and 
increasing with energy and N part , [10]. So we do expect changes in our 
parameters in (23) when moving from central (Au-Au, 0-6%) to peripherical 
(Au-Au, 35-45%) collisions. Our parameterization for peripherical collisions 
gives a' = 0.228 ± 0.002, a = 0.235 ± 0.008 and 5 = 0.90 ± 0.03, to be 
compared with the values in the Table, for central collisions. At much higher 
energy we expect a', a and 5 to approach the central collision values and to 
become the same for all centralities, as well as for pp collisions. 

The parameter Ao = — ay b plays, in our model, the important role of 
controlling the separation (sharply in Models I and II) of the mid rapidity 
dense central region - where percolation dominates - from the fragmentation 
region where, in general, the participating nucleons retain their individual- 
ity. The central rapidity region A central = -ay b - (-y b ) = (1 - a)y b grows 
with the energy. But the fragmentation region, Af ragm . = — (—ay b ) = ay b 
also increases with the energy. This is not the case of an extended plateau 
(Feynmann- Wilson plateau), followed by a fixed rapidity length fragmenta- 
tion region. 

We finally note that the experimental value found for A, not very different 
from the values of Model I and Model II, 0.36, is consistent with values found 
for the intercept of the Pomeron in colour glass saturation models extended 
to AA scattering via geometrical scaling, [16]. 
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Table 



The parameters a', a and S of Eq. (23), in the case of Models I and II, 
and as obtained from PHOBOS data, [2]. 





a' 


a 


5 


Model I 


2/7 


3/7 





Model II 


1/3 


1/3 





RHIC data [2] 








y/s = 19.6 GeV 


0.26 ± 0.002 


0.27 ±0.03 


0.65 ±0.05 


l/s = 62.4 GeV 


0.230 ±0.008 


0.23 ±0.02 


0.62 ±0.07 


y/s = 130 GeV 


0.249 ± 0.006 


0.28 ±0.02 


0.67 ±0.07 


y/s = 200 GeV 


0.251 ±0.005 


0.29 ±0.01 


0.70 ±0.06 


Overall fit to RHIC data 


0.247 ±0.003 


0.269 ±0.007 


0.67 ±0.03 
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